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Introduction 


•  In  an  earlier  paper^,  the  expression  for  the  partial  derivative 
of  an  eigenvalue  with  respect  to  a  parameter  was  given  as 

® k]  ,  (1) 

for  the  eigen-problem  ^  , 

«•„  A  (f„  -  B  ‘P-  =  ^ 

where  A  and  B  need  not  be  Hermitian,  but  distinct  eigenvalues  are  required. 

When  the  flutter  analysis  is  formulated  as  an  eigen-problem,  the 
expression  for  the  derivative  is  applicable.  Both  modal  and  collocation 
flutter  analysis  methods  are  accommodated.  However,  in  the  case  of  modal 
analysis,  it  happens  that  some  of  the  required  quantities  are  not  directly 
available  and  must  be  developed.  The  development  is  relatively  straight¬ 
forward,  and  is  based  on  the  concept  of  completeness  in  an  engineering 
sense.  With  the  partials  of  the  eigenvalue,  the  gradient  of  the  flutter 
stability  constraint  surface  in  parameter  space  may  be  developed,  thereby 
making  possible  automated  optimized  design  which  accommodates  the  flutter 
constraint. 

Modal  Analysis  Method 

In  order  to  provide  easy  reference  and  to  establish  notation,  some 
elementary  results  will  be  given.  The  system  is  represented  by  the  matrix 
equation 


M  X  X  VC  =  9  , 


(3) 


and  the  assumed  homogeneous  solution 

,,  tjnt 

X  =  'If'  e 

leads  to  the  eigen-problem 


Making  the  coordinate  transformation 


'  1 


in  equation  (3)  and  premultiplying  gives 


=  °  ’  (7, 

where  ^  is  a  matrix  whose  columns  are  the  eigenvectors  It  is  noted  that 

the  coordinate  transformation  (6)  assumes  what  is  called  here  completeness 
in  an  engineering  sense.  Under  restriction  of  sinusoidal  motion 


c,-.  (f  si 


m  to 


the  generalized  forces  are  written 


[q]  =  -  [p]  h] , 


where  the  matrix  P  is  a  function  of  reduced  frequency  as  well  as  other 

aerod3mamic  variables. 

The  usual  procedure  is  to  write  the  problem 

■!?]  '^3  +  T’l/P  =  >  (10) 

where  the  g  represents  the  amount  of  structural  damping  required  to  maintain 
sinusoidal  motion.  Equation  (10)  is  in  the  form  of  (2).  In  practice  the 
aerodynamic  coefficients  are  evaluated  at  several  reduced  frequencies,  the 
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corresponding  eigenvalues  found,  and  a  plot  made  of  the  damping  required 
versus  speed.  The  critical  flutter  speed  is  found  by  the  required  damping 
being  zero.  It  is  assumed  that  this  has  been  done,  and  that  the  aerodynamic 
matrix  Pis  evaluated  at  the  critical  speed  or  reduced  frequency.  Obviously 
the' partial  of  the  eigenvalue  of  (10)  is 


u;® 


Proceeding  with  the  development  leads  to 


and  ^ 

-B  =  +  ■P.  . 

yk  V-L- 

The  expression  for  the  derivative  of  an  eigenvector  is  given  by  Fox  and 


Kapoor  as 


where 


'W  ~  T  ^ 


Ej, _ =  ,  (16) 

provided  the  eigenvectors  have  been  normalized  such  that 

'i/'I  M  =  I.  (17) 

\ 

Use  of  equation  (14)  assumes  completeness  in  an  engineering  sense  in 
the  same  way  as  equation  (6).  The  m’"’^  element  of  'ijl  is  .  Therefore,  : 


follows  that 


Therefore,  it 
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or  in  matrix  form 


=  wi>:. 


consider  an  element  P j  of  the  aerodynamic  matrix  F  contained  in 


the  expression 


Zj  '  f  'V-  ‘‘Fj  ’ 


and  take  the  partial  derivative 

2  ■<.  Pj  *  P  =  S •  (21) 

The  pressure  difference  AjOj  caused  by  unit  sinusoidal  amplitude  in  the 
generalized  coordinate  corresponding  to  the  eigenvector  is  found  by 
solving  the  equation  relating  tt  to  the  down  wash  magnitude 

>4  -  ff  ^  ^Pj  »  (22) 

where  K  is  the  known  kernel  function.  Implicit  in  the  development  here  i; 
the  usual  assumption  that  the  pressure  due  to  a  general  downwash  may  be 
written  as  a  linear  combination  of  ^pj .  Taking  the  partial  of  (22)  gives 

^  ‘  ^ ^  ,  (23) 

or,  by  using  (14)  and  (22) 

.‘^,k  *  ^kj«  (24) 

Assuming  an  expression  in  the  convenient  form 

f  (25) 

\ 

(Justified  again  by  assumed  engineering  completeness)  and  writing 


■Z.  =  ff  ^Pj 


then  substituting  into  (23)  yields 


which,  by  using  (22),  reduces  to 

^  Z  V^.j  I  ^  .  (28) 

Equating  (24)  and  (28)  and  simplifying  leads  to 

t 

4  ^  ^  9  ^  4*-^  .  O  .  /0O^ 

from  which  the  may  be  found  by  equating  at  specific  points,  and 

solving  in  matrix  form.  Expression  (21)  becomes 

*^m,k  ^j.k  ”  ^ ^  ( Hja  5 

from  which,  in  matrix  form, 

Ph  ■  **^“*,l<  -  (30) 


[e;-^  . 


Collecting  the  above  results,  equation  (1)  becomes 

r'^V'^3  oj'k'  01/ . 


-  **\v,  k  ^  ’ 


where 


11'?- 

(31) 

■T'. 

(32) 

Numerical  evaluation  of  (31)  may  be  expressed  as 


d^t  i  dj  -t-  UJ„  (t>„  *  '>■  . 


Equating  real  parts  of  (11)  and  (33),  and  solving  gives 

■*  .  «  ' 

“  ~  4.  •<■-43  ’ 

and  similarily,  the  imaginary  parts  (with  g^=  0)  lead  to 

~  (^x  *  4  )  . 
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Collocation  Analysis  Method 


written 


The  eigen— problem  for  the  collocation  flutter  analysis  method  is 
3 


U  =  C[M  vr3k  , 

where  C  is  the  structural  influence  coefficient  matrix,  and!?  is  a 


(36) 


complex  matrix  of  aerodynamic  influence  coefficients,  which  are  functions 
of  reduced  frequency.  Again,  it  is  assumed  that  the  aerodynamic  matrix  is 
evaluated  at  the  critical  condition. 

Proceeding  directly  to  evaluation  of  equation  (1) ,  the  derivative  of 
the  eigenvalue  is 

--  C  k*  ^  .  on 

where 

eZ  ^  ^-1  .  (38) 

It  may  be  easier  to  find  directly  from  the  eigen-problem 

rather  than  to  perform  the  inversion  of  the  large  matrix  in  (38) •  Expressions 
analogous  to  (33),  (34),  and  (35)  maybe  developed. 


Implementation  of  Frequency  Ratio  Limi_t 
It  is  commonly  known  that  flutter  speed  is  sensitive  to  the  ratio  of 
the  two  structural  natural  frequencies  which  couple  to  produce  flutter. 
Johnson  and  Warren^  have  suggested  that  a  limit  be  placed  on  the  ratio  to 
approximate  the  flutter  constraint.  Flutter  analyses  during  preliminary 
design  establish  the  frequency  ratio  limit,  When  the  eigen-problem  of 
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(40) 


the  structure  using  the  structural  influence  coefficient  is  written 

-A  CMOf'  ■<-'1}^  =  o  , 


the  frequency  ratio  limit  approximation  to  the  ]flutter  constraint  may  be 


expires  sed  as 


It  follows  that 


where 


Syk 

2 

and  where  may  be  computed  using  the  Fox  and  Kapoor  expressions  with 

~  ~  ^  ?  (44) 

or  perhaps  the  may  be  taken  equal  to  zero  as  suggested  by  Johnson  and 
Warren.  A  gradient  may  be  developed  from  (42)  and  used  in  an  automated 


design  procedure. 

Discussion 

From  the  physical  meaning  of  the  symbol^,  it  is  clear  that,  at  the 
critical  airspeed,  a  qualitative  improvement  is  indicated  by  the  quantity 
being  negative.  This  indicates  that,  in  effect,  structural  damping 
is  added  to  the  system  by  a  positive  increment  in  the  design  variable  under 
consideration.'  Comparison  of  the  derivatives  of  the  damping  with  respect 
to  the  various  parameters  will  give  some  indication  of  their  relative 
effects.  It  remains  to  place  matters  on  a  quantitive  basis. 
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The  eigenvalue  of  interest  is  a  function  of  the  reduced  frequency  and 
of  the  design  variables,  or  alternatively,  of  the  airspeed  and  design 
variables.  Then  the  increment  in  the  eigenvalue  due  to  increments  in  the 
design  variables  and  in  airspeed  is  approximated  by 

(45) 

To  maintain  the  critical  condition,  the  increment  in  the  imaginary 
part  must  be  zero  (also  0) 


^5-, 


k  ^IPh 


^9-  = 


from  which 


where  the  variation  of  damping  required  as  a  function  of  airspeed  is  obtained 
from  the  flutter  analyses.  Equation  (47)  may  be  used  to  develop  the  gradient 
to  the  constraint  and  to  select  a  design  variable  increment  size. 

In  practice,  it  may  be  more  efficient  to  use  an  approximate  constraint 
(guided  by  previous  flutter  analyses)  such  as  the  frequency  ratio  of 
Johnson  and  Warren^  to  approach  the  optimum  point,  then  to  use  the  present 


method  to  finalize  the  design. 

In  the  development  of  equation  (1),  the  assumption  of  distinct  eigenvalues 
is  made.  Since  there  is  often  a  coalescence  of  frequencies  as  the  critical 
flutter  point  is  approached,  there  may  be  a  practical  numerical  problem 


caused  by  close  eigenvalues. 

Because  of  the  appeal  to  engineering  intuition  in  the  development,  there 
is  a  need  to  verify  the  results  of  this  paper  by  analyzing  typical  problems. 
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Summary 

Expressions  developed  for  the  derivative  of  the  eigenvalue  in  the 
flutter  analysis  problem  offer  potential  both  for  automated  optimum 
design  and  for  less  sophisticated  design  iteration. 
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